***********************Approach 1************************************
 
*sigmaCH -all //signrank sign //
summ sigma0, detail
summ sigma1, detail
signtest sigma0=sigma1
signrank sigma0=sigma1 
 
*sigmaNCH -all // both sign // *REMAINS AFTER EXCUDING large sigmas
summ sigma0n, detail
summ sigma1n, detail
signtest sigma0n=sigma1n
signrank sigma0n=sigma1n
signtest sigma0n=sigma1n if sigma0n<10000&sigma1n<10000
signrank sigma0n=sigma1n if sigma0n<10000&sigma1n<10000

*alphaCH -all
summ alpha0, detail
summ alpha1, detail
signtest alpha0=alpha1
signrank alpha0=alpha1

*alphaNCH -all // signtest sign, signrank 0.07//
summ alpha0n, detail
summ alpha1n, detail
signtest alpha0n=alpha1n
signrank alpha0n=alpha1n

*if subject!=47 (w nchall)
*sigmaCH -all -H
summ sigmah0 , detail
summ sigmah1 , detail
signtest sigmah0=sigmah1 
signrank sigmah0=sigmah1 
 signtest sigmah0=sigmah1 if sigmah0<10000&sigmah1<10000
signrank sigmah0=sigmah1 if sigmah0<10000&sigmah1<10000 

*sigmaNCH -all -H // both sign //// *REMAINS AFTER EXCUDING large sigmas
summ sigmah0n if subject!=47, detail
summ sigmah1n if subject!=47, detail
signtest sigmah0n=sigmah1n if subject!=47
signrank sigmah0n=sigmah1n if subject!=47
signtest sigmah0n=sigmah1n if subject!=47&sigmah0n<10000&sigmah1n<10000
signrank sigmah0n=sigmah1n if subject!=47&sigmah0n<10000&sigmah1n<10000


*sigmaCH -all -T
summ sigmat0 , detail
summ sigmat1 , detail
signtest sigmat0=sigmat1 
signrank sigmat0=sigmat1  
signtest sigmat0=sigmat1 if sigmat0<10000&sigmat1<10000
signrank sigmat0=sigmat1 if sigmat0<10000&sigmat1<10000 

*sigmaNCH -all -T // both sign // *REMAINS AFTER EXCUDING large sigmas
summ sigmat0n if subject!=47, detail
summ sigmat1n if subject!=47, detail
signtest sigmat0n=sigmat1n if subject!=47
signrank sigmat0n=sigmat1n if subject!=47
signtest sigmat0n=sigmat1n if subject!=47&sigmat0n<10000&sigmat1n<10000
signrank sigmat0n=sigmat1n if subject!=47&sigmat0n<10000&sigmat1n<10000


*roCH -all -T
summ ro0 , detail
summ ro1 , detail
signtest ro0=ro1 
signrank ro0=ro1 

*roNCH -all -T
summ ro0n if subject!=47, detail
summ ro1n if subject!=47, detail
signtest ro0n=ro1n if subject!=47
signrank ro0n=ro1n if subject!=47


 


*********************************Approach 2******************************

*sigmaCH -all //signrank sign //
summ sigma0, detail
summ sigma1, detail
signtest sigma0=sigma1
signrank sigma0=sigma1 

*sigmaNCH -all // both sign // REMAINS after excl > 10000
summ sigma0n, detail
summ sigma1n, detail
signtest sigma0n=sigma1n
signrank sigma0n=sigma1n
signtest sigma0n=sigma1n if sigma0n<10000&sigma1n<10000
signrank sigma0n=sigma1n if sigma0n<10000&sigma1n<10000

*alphaCH -all
summ alpha0, detail
summ alpha1, detail
signtest alpha0=alpha1
signrank alpha0=alpha1

*alphaNCH -all // signtest sign, signrank 0.07//
summ alpha0n, detail
summ alpha1n, detail
signtest alpha0n=alpha1n
signrank alpha0n=alpha1n

*if subject!=47 (w nchall and chall)
*sigmaCH -all -H //signtest significant XXXX// -REMAINS after excl > 10000
summ sigmah0 if subject!=47, detail
summ sigmah1 if subject!=47, detail
signtest sigmah0=sigmah1  if subject!=47
signrank sigmah0=sigmah1   if subject!=47
signtest sigmah0=sigmah1 if subject!=47&sigmah0<10000&sigmah1<10000
signrank sigmah0=sigmah1 if subject!=47&sigmah0<10000&sigmah1<10000
 
*sigmaNCH -all -H // signrank sign, signtest 0.07 xxx// - NO result after removing > 10000
summ sigmah0n if subject!=47, detail
summ sigmah1n if subject!=47, detail
signtest sigmah0n=sigmah1n if subject!=47
signrank sigmah0n=sigmah1n if subject!=47
signtest sigmah0n=sigmah1n if subject!=47&sigmah0n<10000&sigmah1n<10000
signrank sigmah0n=sigmah1n if subject!=47&sigmah0n<10000&sigmah1n<10000

*sigmaCH -all -T
summ sigmat0 if subject!=47, detail
summ sigmat1 if subject!=47, detail
signtest sigmat0=sigmat1  if subject!=47
signrank sigmat0=sigmat1  if subject!=47
 
*sigmaNCH -all -T // signrank sign, signtest 0.07 xxx // 0.0557 and 0.19 after removing
summ sigmat0n if subject!=47, detail
summ sigmat1n if subject!=47, detail
signtest sigmat0n=sigmat1n if subject!=47
signrank sigmat0n=sigmat1n if subject!=47
signtest sigmat0n=sigmat1n if subject!=47&sigmat0n<10000&sigmat1n<10000
signrank sigmat0n=sigmat1n if subject!=47&sigmat0n<10000&sigmat1n<10000

*roCH -all -T
summ ro0 if subject!=47, detail
summ ro1 if subject!=47, detail
signtest ro0=ro1 if subject!=47
signrank ro0=ro1 if subject!=47

*roNCH -all -T 
summ ro0n if subject!=47, detail
summ ro1n if subject!=47, detail
signtest ro0n=ro1n if subject!=47
signrank ro0n=ro1n if subject!=47

 
 
  
*********************************Approach 3******************************
  
*sigma-all //signrank sign //
summ sigma0, detail
summ sigma1, detail
signtest sigma0=sigma1
signrank sigma0=sigma1 


*alpha -all
summ alpha0, detail
summ alpha1, detail
signtest alpha0=alpha1
signrank alpha0=alpha1




*sigma -all -H <<signrank sign>> REMAINs after removing > 10000
summ sigmah0 , detail
summ sigmah1 , detail
signtest sigmah0=sigmah1 if sigmah0<10000&sigmah1<10000
signrank sigmah0=sigmah1  if sigmah0<10000&sigmah1<10000


*sigma -all -T
summ sigmat0 , detail
summ sigmat1 , detail
signtest sigmat0=sigmat1 
signrank sigmat0=sigmat1  

*ro -all -T
summ ro0 , detail
summ ro1 , detail
signtest ro0=ro1 
signrank ro0=ro1 

***************************BIASES
*BY CHALL (as in approaches 1 and 2)

*CHALL
 ttest bias0healthier = bias1healthier
 ttest bias0healthier = 0 
  ttest bias1healthier = 0
  ttest bias0rl = bias1rl
  ttest bias0rl = 0
   ttest bias1rl = 0
   
   *biasesHealthier sign. > 0 and p = 0.058 for HL vs LL (two sided)
   
 *NCHALL
  ttest bias0nhealthier = bias1nhealthier
 ttest bias0nhealthier = 0 
  ttest bias1nhealthier = 0
  ttest bias0nrl = bias1nrl
  ttest bias0nrl = 0
   ttest bias1nrl = 0
   *biasesHealthier sign. > 0 in HL only p =0.001
   *biasRight < 0 in LL sign p = 0.045 (two sided)
   
   **No sign differences for comparisons between chall and non-chall
 
*ALL TRIALS (as in approach 3)
 ttest bias0healthier = bias1healthier
 ttest bias0healthier = 0 
  ttest bias1healthier = 0
   ttest bias0tastier = bias1tastier
 ttest bias0tastier = 0 
  ttest bias1tastier = 0
  
  ttest bias0rl = bias1rl
  ttest bias0rl = 0
   ttest bias1rl = 0
   
    *biasRight < 0 in LL sign p = 0.06 (two sided)
	 *biasesHealthier sign. > 0 and p = 0.07 for HL vs LL (two sided)
	 
	 **BOTH 1&2 and 3 (trends): stronger bias to the right (or weaker bias to the left) under load 
	 ***Approach 3 (trend): biases to less tasty